Bethe Ansatz Solution of the Asymmetric Exclusion Process with Open Boundaries 
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We derive the Bethe ansatz equations describing the complete spectrum of the transition matrix 
of the partially asymmetric exclusion process with the most general open boundary conditions. 
For totally asymmetric diffusion we calculate the spectral gap, which characterizes the approach 
to stationarity at large times. We observe boundary induced crossovers in and between massive, 
diffusive and KPZ scaling regimes. 
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The partially asymmetric simple exclusion process 
(PASEP) describes the asymmetric diffusion of particles 
along a one-dimensional chain with L sites. It is one 
of the most studied models of non-equilibrium statisti- 
cal mechanics, see [1, 2] for recent reviews. This is in 
part due to the fact that is one of the simplest lattice gas 
models, but also because of its applicability to molecular 
diffusion in zeolites [3], bioploymers [4], traffic flow [5] 
and other one-dimensional complex systems [6]. 

At large times the PASEP exhibits a relaxation to- 
wards a nonequilibrium stationary state. An interesting 
feature of the PASEP is the presence of boundary induced 
phase transitions [7]. In particular, in an open system 
with two boundaries at which particles are injected and 
extracted with given rates, the bulk behaviour in the sta- 
tionary state is strongly dependent on the injection and 
extraction rates. Over the last decade many stationary 
state properties of the PASEP with open boundaries have 
been determined exactly [1, 2, 8-11]. On the other hand, 
much less is known about its dynamics. This is in con- 
trast to the PASEP on a ring for which exact results 
using Bethe's ansatz have been available for a long time 
[12, 13]. For open boundaries there have been several 
studies of dynamical properties based mainly on numer- 
ical and phenomenological methods [14, 15]. In this Let- 
ter we employ Bethe's ansatz to obtain exact results for 
the approach to stationarity at large times in the PASEP 
with open boundaries. Upon varying the boundary rates, 
we find crossovers in massive regions, with dynamic ex- 
ponents z = 0, and between massive and scaling regions 
with diffusive (z = 2) and KPZ (z = 3/2) behaviour. 

The dynamical rules of the PASEP are as follows. At 
any given time t each site is either occupied by a particle 
or empty and the system evolves subject to the follow- 
ing rules. In the bulk of the system (i = 2, . . . , L — 1) a 
particle attempts to hop one site to the right with rate p 
and one site to the left with rate q. The hop is executed 
unless the neighbouring site is occupied, in which case 
nothing happens. On the first and last sites these rules 
are modified. If site i = 1 is empty, a particle may enter 
the system with rate a. If on the other hand site 1 is 



occupied by a particle, the latter will leave the system 
with rate 7. Similarly, at i = L particles are injected 
and extracted with rates 5 and (3 respectively. With ev- 
ery site i we associate a Boolean variable Tj, indicating 
whether a particle is present (7$ = 1) or not (r, = 0). 
The state of the system at time t is then characterized 
by the probability distribution P t (ri, . . . , tl). The time 
evolution of Pt occurs according to the aforementioned 
rules and as a result is subject to the master equation 



dt 



MP*. 



(1) 



Here M is the PASEP transition matrix whose eigen- 
values have non-positive real parts. The large time be- 
haviour of the PASEP is dominated by the eigenstates 
of M with the largest real parts of the corresponding 
eigenvalues. 

Bethe 's Ansatz: It is well known that the transition 
matrix M is related to the Hamiltonian H of the open 
spin-1/2 XXZ quantum spin chain through a similarity 



transformation M 
by [10] 



pqU\HU x where H is given 



L-l 



h(a z J+1 ~ a*) + A] + + B L . (2) 



The parameters A and h, and the boundary terms -Bi,l 
are related to the PASEP transition rates by 




Br = 



P + S-tf- 5)aj - - 25XQ L - 1 al 

2VM 



a + 7 + {a - i)a\ - 2aAaf - ^cr+ 



(3) 



Here A is a free parameter on which the spectrum does 
not depend and cr^ = (er^ ± icrJ)/2. 
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Although it has been known for a long time that H is 
integrable [16], the off-diagonal terms in B\ and Bl have 
presented great difficulties in diagonalizing H using e.g. 
Bcthe's ansatz. However, recently a breakthrough was 
achieved [17] in the case where the parameters satisfy a 
constraint, which in our notation reads 



(Q L+2k - l)(a/3 - jSQ 



L-2fc-2\ 



= 0. 



(4) 



Here k is an integer such that \k\ < L/2. For a given k 
this constraint can be satisfied by choosing Q to be an 
appropriate root of unity, or by relating the boundary 
and bulk parameters such that the second factor in (4) is 
zero. However, for generic values of the PASEP parame- 
ters (4) can also be satisfied by choosing k = —L/2. For 
this choice of k we infer from [18] that for even L there 
is an isolated level with energy Eq = 0, the ground state 
energy of the PASEP. Furthermore, all excited levels are 
given by 



L-l 



E = a + l3 + 1 + 5+Y, 



(Q 2 



3 = 1 



(Q- Zj )(Q Zj -iy 



(5) 



where the complex numbers Zj satisfy the Bethe ansatz 
equations 



ZjQ 



Q 



2L 



L-l 



K( Zj ) = n 



Zl ZjZiL 



1 



Zj - ziQ 2 zjzi - Q 2 



Here K(z) = K(z, a, j)K(z, (3, 5) and 

-az 2 + Qz(Q 2 - 1 + a - 7) + 7Q 2 



~/Q 2 z 2 + Qz(Q 2 - 1 + a - 7) - a 



(6) 



(7) 



In order to ease notations we have set, without loss of 
generality, p = 1 and hence Q = ^/q. We note that 
in the case of symmetric diffusion Q = 1 (6) reduce to 
the Bethe ansatz equations derived in [19] by completely 
different means. In order to determine the exact value 
of the spectral gap we have analyzed (5) and (6) in the 
limit L — ► 00. To simplify the analysis, we will focus on 
the case of total asymmetry 7 = <5 = 0,Q^0in the 
remainder of this Letter. 

Totally asymmetric exclusion (TASEP): After a rescal- 
ing z — > Qz and setting 7 = S = 0, the Q — > limit of 
equations (5) and (6) reads 



L-l 



E 



+ / ? + E^ L T 

z l - 1 



1=1 



( Zj If 



L-l 

{ Zj + a) {zj + b) ]J (zj 

1^3 



(8) 



, (9) 



where a = (1 — a) /a and b = (1 — /?)//? . We define 
g(z) = mz/(z — l) 2 and t/b (z) = lnz/(l — z 2 )+ln (z + a) + 
In (z + b) , and consider the "counting function" [20] , 



1 1 

\Y L (z) = g(z) + -g h (z) + - ^2 K(z u z), 

1=1 



(10) 



where K(w, z) — — mw + ln(l — wz). Equations (9) can 
now be written as 

Y L { Zj ) = T I 3 C ? = 1,...,L-1), (11) 

where the Ij are integers. Each set of integers {Ij} spec- 
ifies a particular excited state, and we find numerically 
that the first excited state is obtained for the choice 

Ij = -L/2 + j for j = l,..., L-l. (12) 

The eigenvalue (8) can be written as 

E = a + (3 + L\im(iYi{z)-g'{z)-y h {z)} . (13) 

In order to derive an integral equation for Yl(z) in the 
limit L — > 00 we write, 

(14) 

where C = C\ + C2 is a contour enclosing all the roots Zj , 
Ci being the "interior" and Ci the "exterior" part, see 
Fig. 1 . The contours C\ and C2 intersect in appropriately 
chosen points £ and £*. Using the fact that integration 




FIG. 1: Sketch of the contour of integration C in (14). The 
open dots correspond to the roots Zj and £ is chosen close to 
zl-i and avoiding poles of cot(LYt(x)/2). 

from £* to £ over the contour formed by the roots is equal 
to half that over C2 — C\ we find using (14), 

iY L (z) = g(z) + j 9h (z) + i- J^K(w,z)Yl(w)dw 



1 r k 



K(w,z)Y[(w) 



iLY L (w) 



+ 



K(w,z)Y[(w) 



2tt Jc 2 e iLY ^( w )-l 



dw 



dw, 



(15) 



where we have chosen the branch cut of K(w, z) to lie 
along the negative real axis. The eigenvalue (13) is ob- 
tained by systematically expanding (15) in the system 
size using standard methods. We note that care has to 
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be taken when there is a stationary point of Yl(z) close 
to the contour of integration, in which case an analysis 
similar to that in [13] has to be carried out. 

Let us briefly recall the stationary state phase diagram 
derived in [8, 9]. There are altogether four phases in the 
stationary state at t = oo: (1) the low density phase 
for a < (3 < 1/2; (2) the high density phase for j3 < 
a < 1/2; (3) the coexistence line at (3 = a < 1/2; (4) the 
maximal current phase at a, (3 > 1/2. We now determine 
the scaling of the spectral gap in these regimes from the 
Bethe ansatz equations. 

Low and High Density Phases: Let us fix the end points 
r and £ by 



TT 



y l {0 



TT 

T 



(16) 



The integral over C\ in (15) can be calculated by splitting 
the contour into its upper and lower parts and expanding 
the integrand around £ and £* respectively. Expanding 
in inverse powers of L, i.e., 

Y L ( Z ) = J2 L ~ nY n(z), Z = Z c + J2*nL- n , (17) 



n=0 



71=1 



and assuming that — a < z c and —b < z c , we find from 
equation (15) to C(L _1 ) that, 



Yo(z) = 
Yi(z) = 



-iln 



- i In 



z 
z c 

z 1 



1 - Zc 

1 - z 

-z 2 



(18) 



Z c 1 - Z 2 

-iln {ab(-ZcD, 



Z - Z c 



z + a z + b 



Zc + azc 



^b 
(19) 



where v\ = —Yq(z c )5i/tt. The values of v\ and z c follow 
from (16) to be v\ = 2, z c — —l/y/ab. Substituting these 
values into (13) we obtain the gap (20), which is of order 
0(1) in the limit L — > oo. 

If —b > —l/Vab the point —6 lies inside the contour 
formed by the roots, see Fig 1, giving rise to a different 
solution for Y\(z). Comparing again with condition (16) 
we find in this case that v\ = 3 and z c = — a -1 / 3 = 

— l/^abc, resulting in the spectral gap given in (21), 
which is independent of (3 and again of order 0(1) in 
the limit L — > oo. A similar analysis is made when 
—a > —l/\/ab. As the spectral gap is 0(1) in the low 
and high density phases, the correlation length is finite 
and these phases are therefore massive. 

Coexistence Line: Subleading corrections can be ob- 
tained by taking higher order terms into account in (17). 
After a lengthy calculation we find that the gap vanishes 
like L~ 2 on the coexistence line (3 = a, with a constant 
of proportionality given by (22). 

Maximal Current phase: When (3 > (3 C and a — ► 1/2, 
the value of z c where the contour closes approaches z c — 

— 1 and a straightforward expansion of the last two terms 
in (15) breaks down as Y[(C) ~ Y^(-l) = 0. A further 



complication is the singularity in K(w, z) at w = z = 
z c = — 1. To proceed one expands around z c defined by 
Y[( z c) = [13, 21]. This gives rise to an expansion of 
Yl{z) in powers of L' 1 / 2 and one finds in lowest order 
that the energy gap vanishes as L^l 2 . The prefactor can 
only be determined numerically. 

We now summarize our results. We have used Bethc's 
ansatz to diagonalize the PASEP transition matrix M for 
arbitrary values of the rates p, q, a, (3, 7 and S that char- 
acterize the most general PASEP with open boundaries. 
The resulting Bethe ansatz equations (5), (6) describe 
the complete excitation spectrum of M and are one of 
our main results. We have carried out a detailed analy- 
sis of the Bethe ansatz equations for the simplified case 
of the TASEP and determined the exact asymptotic be- 
haviour of the spectral gap for large lattice lengths L. 
This gap determines the long time (t ^> L) dynamical 
behaviour of the TASEP. We emphasize that care has to 
be taken regarding time scales, and that our results be- 
low are not valid at intermediate times t « L where the 
system behaves as for periodic boundary conditions [2]. 

We found that there are three regions in parameter 
space where the spectral gap is finite and the station- 
ary state is approached exponentially fast, and one re- 
gion and a line where the gap vanishes as I — ► 00. 
The resulting dynamical phase diagram is shown in Fig- 
ure 2. In order to parametrize the phases we define 
(3 C = (1 + a- 1 / 3 )- 1 and a c = (1 + 6" 1 / 3 )- 1 . The values 
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FIG. 2: Dynamic phase diagram of the TASEP. Mi, M na and 
Mub are massive phases, CL denotes the critical coexistence 
line and MC the critical maximal current phase. 

of the spectral gap in the various regions of the phase 
diagram of the TASEP are as follows: 
Massive Phase Mj: a < a c , (3 < (3 C , a 7^ (3 

-E 1 =a + P-—l T =+0(L-% (20) 
1 + vab 

The spectral gap does not vanish as L — > 00 and hence 
implies a finite correlation length and exponential ap- 
proach to stationarity. 
Low Density Phase Mn a : a < 1/2, (3 > (3 C 

2 



-E 1= a + I3 C 



(21) 
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Note that in this phase the spectral gap is independent 
of 0. The behaviour in the high-density phase Mnt>: < 
1/2, a c < a is obtained by the exchange a <-> 0. 
Coexistence Line (CL): = a < 1/2. 

-E 1 = ^-^L^ + 0(L-% (22) 
1 — 2a 

We thus find a dynamic exponent z = 2 corresponding 

to diffusive behaviour. 

Maximal Current Phase (MC): a,0> 1/2. 

-E x S3 3.56L" 3 / 2 + 0(L- 2 ). (23) 

In this phase, which coincides with the stationary max- 
imum current phase, we find a KPZ [22] dynamic expo- 
nent z = 3/2. The gap is smaller than that of the periodic 
case where it is found that -E 1 = 6.509 . . . L^/ 2 [13, 23]. 
We note that the subdivision of the massive high and low 
density phases is different from the one suggested on the 
basis of stationary state properties in [9]. 

Discussion: It is known [24] that by varying the bulk 
hopping rates one can induce a crossover between a dif- 
fusive Edwards- Wilkinson (EW) scaling regime [25] with 
dynamic exponent z = 2 and a KPZ regime [22] with 
z = 3/2. In this letter we have shown using exact meth- 
ods that a crossover between phases with z — 2 and 
z = 3/2 occurs in the case where the bulk transition rates 
are kept constant, but the boundary injection/extraction 
rates are varied. As shown in [15] the diffusive relax- 
ation (z = 2) is of a different nature than in the EW 
regime and is in fact due to the unbiased random walk 
behaviour of a shock (domain wall between a low and 
high density region). Our results (20) and (22) for the 
massive phase Mi and the coexistence line agree with the 
relaxation time calculated in the framework of a domain 
wall theory (DWT) model in [15]. This is in contrast 
to the massive phases Mn, where the exact result (21) 
differs from the DWT prediction. An interesting open 
question is whether it is possible to understand (21) in a 
generalized DWT framework. 

The Bethe ansatz equations (5), (6) allow for the exact 
determination of further spectral gaps. We find that the 
eigenvalue of the transition matrix with the next largest 
real part is complex, which leads to interesting oscillatory 
behaviour at large times. The dynamical phase diagram 
for the general PASEP is expected to be significantly 
richer, and its analysis is under way. 

The condition (4) is a reflection of the non-semisimplic- 
ity of an underlying Temperley-Lieb algebra with two 
additional boundary generators [26]. Remarkably the 
PASEP satisfies this constraint for arbitrary values of its 
parameters. Generically, non-semisimplicity implies cer- 
tain symmetries in the spectrum, and the physical con- 
sequences of these are currently under investigation. 
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